Accepted for publication in PASP 

Preprint typeset using 1^1^^ style cmulatcapj v. 2/19/04 



00 

O 
O 
(N 

< 



6 



> 

(N 

^' 

O 
oo 
O 



- 

X 



CLOSURE PHASE SIGNATURES OF PLANET TRANSIT EVENTS 

Gerard T. van Belle^ 

European Southern Observatory, Garching, Germany 85748 
gerard.van.belle@eso.org 
Accepted for publication in PASP 

ABSTRACT 

Planet transit events present as attractive targets for the ultra-high-resolution capabilities afforded 
by optical interferometers. Herein is presented an evaluation of the possibility of detection of such 
events through measurement of high-precision closure phases with the MIRC instrument on the 
CHARA Array. Recovery of the transit position angle upon the sky appears readily achievable with 
the existing capabilities of the instrument, along with characterization of other system parameters, 
such as stellar radius, planet radius, and other parameters of the transit event. This technique is 
the only one presently available that can provide a transiting planet's orbital plane position angle, 
and can directly determine the planet's radius independent of any outside observations, appearing 
able to improve substantially upon other determinations of that radius. Additional directly observed 
parameters - also not dependent upon transit photometry or spectroscopy - include impact parameter, 
transit ingress time, transit velocity and stellar radius. 

Subject headings: techniques: interferometric, techniques: high angular resolution, stars: planetary 
systems, stars: individual: HD189733 



1. INTRODUCTION 

Recent discoveries of stars exhibiting the telltale signs 
of planet transits has begun to add a new layer of un- 
derstanding to the rapidly developing field of exoplane- 
tology. While the technique of radial velocity detection 
has prod uced the g reatest yield to date of planet de- 
tections ((Butler et al.,.2006.) . the detected transit events 
have served to define the specific nature of those planets, 
including parameters such as density, at mospheric com- 
posit ion, and aspects of system dynamics (iBurrows et al.l 

Advances in the state of the art in astronomical optical 
interferometry can be directed at these recent transit dis- 
coveries and also contribute to filling in the pieces of the 
exoplanet puzzle. Specifically, measurements of interfer- 
ometric closure phase during a planet transit event can 
determine the inclination and orientation of the plane- 
tary orbit upon the sky, in addition to refining the an- 
gular diameter measurements of both the plane t and the 
star. Just as the Ross iter-McLaughlin effect ()Rossitei| 
119241 : lMcLaughlinlll92^ in radial velocity measurements 
can contribute t o our knowledge of transiting planet sys- 
tem parameters (|Winn et al.ll2006f ) , transit event closure 
phases can further the physical description of these sys- 
tems through direct detection of the transit. Interfero- 
metric phase in general is a powerful tool that is begin- 
ning to be expl oited to its fullest potential in astronomy 
(|Monniei]l2f)fm . 

Planet transit closure phase observations described 
herein are the only presently available technique that 
provide a measurement of the transiting planet's orbital 
plane orientation upon the sky. These closure phase ob- 
servations also uniquely determine the other observables 
of the system - impact parameter, transit velocity, stel- 
lar radius, planet radius, transit ingress time - without 
the need for supporting observations such as transit pho- 

^ For preprints, please email gerard.van.belle@eso.org. 



tometry. For example, the previous direct det ermination 
of HD189733b's diameter (jBaines et al.ll2007l ) measured 
that parameter through a combination of interferometric 
measurements and transit photometry; this technique is 
independent of such outside measurements. 

We will begin with a review of instrument capabilities 
in [21 examination of known extrasolar planet candi- 
dates in fj3l 8- 'quick-and-dirty' partial analytic solution 
of the problem in 21 3' more thorough discussion of clo- 
sure phase leading to a numeric model and analysis of 
planet transits in ^ and finally a full Monte Carlo sim- 
ulation to recover synthetic transits in |j6l 

2. INSTRUMENT CAPABILITIES 

The Georgia State University Center for High An- 
gular Resolution Astronomy (CHARA) Array is a six- 
element optical interferometer located atop Mount Wil- 
son in southern California. The CHARA Array consists 
of six 1-m telescopes laid out on a 'Y' array, two tele- 
scopes per arm, with a baselines of > 300 meters on the 
3 longest baselines. Initial scien ce operations and a fa - 
cilit y description can be found inlMcAlister et al.l ()2005D 
and iten Brummelaar et al.l (|2005l )^ 

Commissioning on-sky tests have recently begun at the 
CHARA Array with the Michigan Infrared Combiner 
MIRC) , a high-precis ion multi-telescope beam combiner 
Monnier et"alll2006af) . MIRC's capability to combine 4 
or 6 telescopes simultaneously and provide 3 to 10 closure 
phase measurements on sources represents a major step 
forward in capability for the facility. Moreover, initial 
MIRC tests are indicating a remarkable ability to mea- 
sure closure phases with precision unprecedented in the 
field of optical interfer ometry, at a level of (T$ ~ 0.03° 
iMonnier et all (|2006af ): the first science demonstrated 
by this capability includes direct imagin g the surface o f 
the rapidly rotating star Altair (Monn ier et al.|[2007f ). 
Other instruments, such as the VLTI AMBER instru- 
ment (jRantakvro et al.ll2006f ). also provide the capability 
to make closure phase measurements, although initial in- 
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dications are that the closure phase precision of AMBER 
is not quite as capa ble as MIRC, with ag, ^few degrees 
(|Weigelt et al.l[200l . 

3. POTENTIAL PLANET TRANSIT TARGETS 

The obvious candidate for o bservations of a 

planet transit even t is HD189733 (jBakos et all l2006t 
iBouchv et al."2005'). With an angular size of 376 ± 31 
/Lias (Baincs et al. 2007), it is the planet-transit hosting 
star with the largest a ngular size dis c overe d to date. 
The discovery paper of iBouchv et all (|2005D cites the 
following system parameters: (a) a planet-star radius 
ratio of 0.172 ± 0.003, (b) an orbital inclination of 
i = 85.79 ± 0.24, (c) an orbital radius of 0.0313 ± 0.0004 
AU, and (d) a transit duration of roughly 1.7 hours. 
The geometry of the transit is depicted in Figure [T] The 
best known values for the host star and planet are found 
in Tables 1 and 2, respectively. 

The next best candidate known at the time of th is draft 
is GJ 436 (|Butler et al.l[200l iGillon et al.l[200l . with 
roughly the same anticipated angular size as HD189733. 
HD149026 (Sato et al.l]2Q 0l, HD17156 (iFischer et al] 
2007"; 'Barbieri et al."2007^, and HD209458 (jHenrv et all 



1999; Charbonncau ct al. 2003) are also worth consid- 
ering, although their stellar angular diameters of 170- 
250 /las are significantly less favorable for detection 
when considering the currently available capabilities of 
CHARA-MIRC. 

4. PARTIAL ANALYTIC SOLUTION FOR TRANSIT 
INTERFEROMETRIC VISIBILITY 

The complex interferomctric visibility of a binary star 
can be written as: 



T/ I ^T/ ^~2'jTi(uAa+vA[3) 



where r is the brightness ratio ()Herbison-Evans et al.l 
[T97lh . Va and Vb are the visibility functions asso- 
ciated with a uniform disk, V — 2Ji(x)/(x), where 
X — Ojjjj'kB I \^ duo is the uniform disk angular diameter, 
B is the projected baseline, and A is the wavelength of 
operation. Dropping the absolute phase term and rewrit- 
ing this in terms of relative separation vector As and 
baseline vector B: 



binary 
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(2) 



Determining visibilities for the specific case of a planet 
transit can be adopted from this formalization with the 
following caveats: 

• The 'brightness ratio' r is the negative value of the 
squared planet-to-star diameter ratio. 

• Equation [2] assumes that both the planet and stel- 
lar disks are indeed uniform disks, since r is con- 
stant. Specifically, limb darkening of the star is 
ignored. 

• Equation [5] is valid only for the portion of the tran- 
sit event when the planet is fully in front of the 
star, when r is again constant. The ingress and 
egress portions of the transit event are not prop- 
erly represented by this equation. 



This solution is useful in characterizing the order- 
of-magnitude effects for observation planning; however, 
given the above caveats, it is insufficient for proper eval- 
uation of actual data. 

5. TRANSIT CLOSURE PHASE AND VISIBILITY 
DEVIATIONS 

In its simplest realization, an optical interferometer 
measures the Fourier components of an image upon the 
sky. The location of the image's components in the 
Fourier transformed {m, v} plane are dictated by the 
baseline between the telescopes in the interferometer pro- 
jected towards the source of interest, and the wavelength 
of operation. As seen in ijH at least a partial analytic 
solution can be predicted from the parameters of the 
experiment. In practice, only the visibility amplitude 
(typically just referred to as the visibility) can be mea- 
sured, while atmospheric turbulence corrupts the direct 
measurement of visibility phase. 

However, interferometers using three or more tele- 
scopes can produce a measurement of the closure phase, 
a phase quantity tha t remains uncorrupted by tele scope- 
specific phase errors (jJenniso onnie l2000h . The 
closure phase $ is the sum of visibility phases around a 
closed loop of baselines. For three telescopes i,j,k, this 
is easily deduced from the observed visibilities: 



^ijk = arg{Vi.j) + arg{V.jk) + arg{Vki) 



(3) 



Each pair of telescopes produces a source visibility and 
phase, where the phase is associated with the source's 
intrinsic phase 0, phase errors 9i — 02 associated with the 
telescope pair, and noise. Typically due to atmospheric 
corruption, for a pair of telescopes, phase information is 
useless. However, for a three telescope array {l,m,n}, 
combination of the three measured phase pairs {ipim — 
<t>im + {Qi — dm), ■ ■ ■) results in cancelation of the phase 
error terms 9 leaving only the sum of the three source 
phases intrinsic to the object - the closure phase, ^imn = 

4>lm + 0mn + 4>nl- 

Use of the closure phase effectively cancels many of 
the corrupting effects of the atmosphere and the instru- 
ment, and is a highly sensitive probe for interferomctric 
image construction on the smallest spatial scales. Signif- 
icantly more complete discussions of the topi c of cl osure 
phase may be found in ' Pearson fc ReadheadI ()1984f ) and 
Monnicr (2007). Closure phases have been used to ex - 
plore disk asymmetries in YSOs (|Monnier "erani2006br ) 
and are very sensitive to asymmetries in images, which 
will prove quite useful in the application discussed here. 

For a star with a planet blocking out part of its disk 
during a transit event, the degree of asymmetry is ex- 
treme - significantly much more so than a star with a 
starspot on its surface: the spot temperature is merely 
some slight fraction of the rest of the photosphere and is 
still emitting radiation at a fairly significant level. At the 
near-infared wavelengths being considered here, a tran- 
siting planet emits extremely little radiation with regards 
to the area of the stellar photosphere it is blocking off 
from our line of sight. 

For the full envelope of expected visibility amplitudes 
and closure phases for a gedanken experiment covering 
the interferometer response during a planet transit event, 
the analytic solution of SjHis insufficient, in that it breaks 
down during the transit ingress and egress. As a means to 
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Fig. 1. — On-sky geometry for a model representative of HD189733, based on data found in lBouchv et al.l 1)20051 ') and IBaines eTan pOOT'l. 
as discussed in JS] In particular, the planet radius of -Rpianot = 64/ias and impact parameter b = 120/ias are derived from the stellar radius, 
the planet-star radius ratio, and orbital inclination. The values of ci = 181 and C2 = 95/ias mark the distance of the beginning and end of 
the ingress event from the transit meridian. The on-sky orientation of the transit chord is, at present, unknown. 



TABLE 1 

Summary of parent star parameters for HD189733. 



Parameter 


Value 


Units 


UD Angular Size 


366 ± 31 


/J as 


Limb darkening coefficient 


1.9 


% 


LD Angular Size 


372 ± 31 


fias 


Bolometric Flux 


28.28 ± 0.49 


10~^ erg cm~^ 


Parallax 


51.94 ±0.87 


mas 


Effective Temperature 


4980 ± 200 


K 


Linear radius 


0.779 ± 0.066 


Rq 


Mass 


0.82 ± 0.05 


Mq 


Luminosity 


0.305 ± 0.008 


Lq 


log g 


4.607 ± 0.043 


[cm s ^] 



Reference 



Baines et al. (2007) 
Tango & Davis ( 20021 
Baines et al. (2007) 
Baines et al. (2007) 
Ferryman et al. (1997) 



Baines et al. (20071 
Bouchy et al. (2005) 



explore that full envelope, a numeric analysis can be per- 
formed to compute directly the expected visibility com- 
ponents. 

For the specific case of HD189733, such an ex- 
periment may easily be executed, leading to an 
expectation of the changes in observed visibility am- 
plitudes and closure phases during the planet transit 
event. For the gedanken experiment, three baselines 



were postulated, with coordinates (in meters) of 
{125.35, 305.94}, {-300.42, -89.62}, {175.07, -216.32} 
(with zero vertical separation), which correspond 
roughly to the CHARA ElSl, WlEl, and SlWl 
baselines, respectively. 

From the sy s tem p arameters for HD189733 as cited by 
iBouchv et al.l ()2005[ ). a model of the transit event was 
constructed with a stellar radius of 185/xas, a planet ra- 
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TABLE 2 

Summary of exoplanet parameters for. HD189733. 



Parameter Value Units Reference 



Radius ratio 0.172 ± 0.003 Bouchy et al. 2005 

Angular diameter 64.0 ± 5.4 fj,as 

Linear Radius 1.19 ± 0.10 Rjup Baines et al. (2007 ) 

Density 0.91 ± 0.23 g cm-^ Baines et al. (2007 ) 



case, but by only a marginal amount - on the order of 
±0.01%. Such a measurement is beyond the capabilities 
of any existing interferometer by two orders of magni- 
tude. Howe ver, the closure phase e xcursion is ±0.2°. As 
detailed in iMonnier et al.l (|2006al ). closure phase mea- 
surements at this level of precision appear possible: ini- 
tial tests of the CHARA-MIRC system showed closure 
phase formal error at the (T$ ~ 0.03° level over the course 
of 3 hours. Shorter integration times indicated a corre- 
spondingly higher level of scatter, but the magnitude of 
this error gives a starting point from which to evaluate 
the possibility for observation of a planet transit event 
using closure phases. 

6. RECOVERY OF TRANSIT PARAMETERS FROM 
CLOSURE PHASE MEASUREMENTS 

Having demonstrated in Sj5]in at least a qualitative way 
that closure phase excursions result from a planet transit 
event, it is useful to further demonstrate that an event 
can be reconstructed from an ensemble of closure phase 
measurements taken during a transit. Free parameters 
of the fit are the primary descriptors of the image upon 
the sky: 

• The stellar radius Tstar- 



dius of 32^as, and a transit chord that was offset from 
the center of the stellar disk by 121/j,as. The rate of the 
transit was not considered in this section, in this simple 
inspection of the effects upon the interferometer visibil- 
ity signals, although in fJH] variables for that aspect of 
the event will be introduced. The top panels of Figure 
[2] shows the this transit event for three different orienta- 
tions upon the sky, a = 0°, 45°, 90°. 

For each of the three orientations, the visibility for each 
of the three baselines, along with the closure phase, was 
computed. This computation was performed in the fol- 
lowing manner: An appropriately limb darkened model 
star was created numerically on a 1024 x 1024 grid with a 
diameter correspondi ng to t he HD189733 parent star. As 
introduced by iMilni (|1921h and discussed in the contex t 
of stellar interferometrv bv lHanburv Brown et al.l (|1974f ) , 
the conventional linear representation of limb darkening 
across the disk of a star can be written as: 

/a(a')-/a(1)[1-«a(1-m)] (4) 

where /i = cos 7, 7 is the angle between the line of sight 
and the stellar surface normal, and /a(1) is the specific 
intensity at t he center of the disk. For HD189733 at 
5000±100K ( Baines et ani2007[ ). the hmb darken ing pa- 
rame ter w at 1.6 /xm is roughly equal to 0.35 (Clare t et al.l 
fl995[) : as we shall see in [JH this technique is relative in- 
sensitive to limb darkening of the planet host star. This 
latter fact is unsurprising given the small angular size 
of the star relative to our notional array; for an array 
with larger baselines {B ^ 1km), there would a greater 
sensitivity to this parameter. 

A fully darkened spot corresponding to the radius of 
the transiting planet was then created on the image for 
a given location along that transit, and the Fourier com- 
ponents were computed. Rather than bear the full com- 
putational load of Fourier transformi ng the entire ima ge 
upon the sky, the approach of .Aufdenberg et al.l (|2006[ ) is 
followed, and only the specific components corresponding 
to the 3 baselines in question were computed, resulting 
in a much lighter computational load without a sacrifice 
in precision. 

This process was then repeated for various points along 
the planet transit. The computed values for visibility 
and closure phase were compared to the nominal values 
for the uneclipsed parent star, and those deviations are 
plotted in the middle and lower panel of Figure [2] This 
approach to computing the visibility amplitudes and clo- 
sure phases can be seen as being superior to the analytic 
solution in [21 since it takes into account stellar limb 
darkening, and is valid through the transit event, includ- 
ing ingress and egress. 

As seen in Figure [21 during the transit event, the vis- 
ibility deviates from from the nominal unocculted star 



• The planet to star radius ratio R. 

• The orbit orientation upon the sky a, defined as 
the angular orientation of the chord of the transit 
event across the disk of the star in right ascension 
and declination, as measured from north to east on 
the sky. 

• The 'impact parameter' 6, defined as the distance 
between the chord of the transit event across the 
disk of the star, and the center of the stellar disk. 

• The 'zero time' JDq of the transit event, defined 
here as the time of closest approach of the planet 
disk to the center of the stellar disk. 

• The velocity v of the planet disk across the disk of 
the star. 

Fixed input parameters of the fit are: 

• The i = 1 ... TV closure phases $i and their as- 
sociated errors cr$_i. For this exercise, it will be 
assumed that these are the single closure phases 
associated with the non-diurnally evolved CHARA 
baselines associated with the SI, S2, El, and 
Wl stations. In practice, CHARA-MIRC is op- 
erated with either four or six simultaneous base- 
lines, which in the latter case can potentially pro- 
vide even more closure phases than being modelled 
here; this, in conjunction with baseline diurnal evo- 
lution, would serve to further constrain the transit 
event parameters. 

• The closure phase errors a^^i. These errors were 
assumed to have a normal distribution, and for 
the various synthetic data sets created below, as- 
sumed to be of a magnitude ranging from cr$ ^ = 
0.005°, 0.010°, 0.020°, 0.050°, 0.100°. The FWHM 
of the cr$.i distribution was set at one-quarter of 
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the average error value. Individual closure phases 
$i were randomized by an error of the scope of CT$,i 
(but independent thereof). 

• The operational wavelength of each observation 
of the closure phases. A single wavelength of 1.6 
/ini will be assumed here for the synthetic data sets. 
In practice, the MIRC instrument spectrally dis- 
perses the starlight and multiple closure phases per 
observation are available as a result. 

• The time ti of observation of each closure phase. 
For our synthetic data sets, we will postulate that 
sets of closure phases are taken sequentially at in- 
tervals of 0.001 of a Julian Day (roughly 85 seconds 
set-to-set), for a total duration of 0.124 days (about 
3 hours), which spans the duration of the transit 
event, plus about 40 minutes before and after the 
transit event. This resulted in iV = 124 data points 
per observation set. 

• The {m, v} coordinates of observing baselines for 
each closure phase, For an actual observation, the 
baselines would need to be properly projected onto 
the sky incorporating diurnal motion. 

• The limb darkening of the star. 

The stellar limb darkening is also potentially fittable free 
parameter; however, for the marginal resolution case of 



CHARA observing HD189733, a model value is sufficient 
for the fitting. 

For simplicity of this investigation, diurnal motion will 
be ignored and the {m, v} coordinates of the observing 
baselines will be constant. In practice, this motion will 
need to be accounted for but will actually provide addi- 
tion constraints upon the image reconstruction, much in 
the same way that baseline evolution can serve to assist 
in constraining the parametrization of binary star orbits 
(|Boden et al.lll999l ). 

To create synthetic 'observation sets' for testing our 
fitting and parameter recovery routines, the free pa- 
rameters {rstar,R,b,JDo,v} discussed at the begin- 
ning of this section were set to the values seen in our 
HD189733-like system in Figure [1] (with JDq = -6000 
sec, R — 0.17, b = 121.0/ias, rgtar = 185^as, and v ~ 
0.07000/ias/sec), and data sets were generated within the 
context of the fixed input parameters discussed above. 
For those data sets, the position angle a of the planet 
transit across the stellar disk was also set to the values 
28°, 115°, and 170°, to test the sensitivity of the param- 
eter recovery on that particular parameter as well. 

Thus, to test 'goodness of fit' for a given set of 
six randomized free parameters {rgtar, R, a, b, JDq, v}, a 
model transit event sequence was generated, projected 
upon the sky, and resultant image sequence Fourier 
transformed for comparison to each of the observed 
(^i data points, and a x^/DOF calculated. A multi- 
dimensional optimization code was then utilized to loop 
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about this goodness of fit routine and derive the best 
{rstar, R,a,b, JDo,v} solution from any given start- 
ing point, a proc ess that took typically 500 iterations 
(jPress et al.lll99^ . An exhaustive search of the transit 
event parameter space was used to explore the x^/DOF 
space, using a two- fold approach. First, a grid of reason- 
able {rstar, R, ct,b, JDo,v} starting values was explored 
to see if the original {rgtar, R,o:,b, JDq,v} parameters 
could be recovered, with a varying over its full range 
of {0°,360°}, and the other parameters being explored 
over a range of ±50% of their 'true' values. These lat- 
ter ranges were expected to encompass the reasonable 
starting points for an actual investigation, based upon 
constraints that may be available from discovery pho- 
tometry or spectroscopy. Each range was gridded with 
a density of 5 to 10 points per variable. Second, a 
large number of iterations {N ~ 1000) were also run 
for each synthetic data set starting from fully random- 
ized {vstar, R, ca, b, JDq, v} Starting values, also with the 
purpose of recovering the original {rgtar, R, ol, 5, JDoi v\ 
parameters that described the generating values of the 
synthetic data set. The resulting manifold appeared 
to be smoothly varying, as the recovery of the original 
parameters appeared to occur without local minima ob- 
structing the recovery of the global minima. Once the 
best solution was established for a given data set, 1 — cr 
errors were established about the minimum through 
exploration of appropriate Ax^ intervals. The results for 
each of the three transit position angles, with the five 
different levels of closure phase error, are seen in Table 

m 

7. DISCUSSION AND CONCLUSION 

Examination of Table [3] illustrates that, even with a 
crude level of closure phase error (ct$ ~ 0.1°), the po- 
sition angle of the transit event is readily recovered to 
within a few degrees, and markedly better with modest 
improvements in closure phase error. Since closure phase 
is, in essence, an observable that quantifies the degree of 
asymmetry in an image upon the sky, we expect that 
this technique should work well even for grazing transit 
events. These events would be limited presumably by the 
shorter duration of the occultation event (and thus fewer 
closure phase data points to fit), and also by a lesser stel- 
lar surface area occulted by the planet's disk, resulting in 
a smaller closure phase signal, but the basic promise of 
the approach holds true. The sensitivity of this approach 
to limb darkening - mentioned, but largely dismissed in 
fJS]- is only slight, due to the largely axisymmetric nature 
of limb darkening. 

Of particular interest is the fact that, in the best 
conceivable case for each apparent transit position an- 
gle, the planet radius appears to be recoverable to the 
level of roughly one part in 30-40, which appears to 
best the previous intcrferomctric measure of the plan- 
etary diameter {Baines ct al. 2007) by a factor of 2.5. 
Additionally, the technique establishes that diameter in 
a w ay that is in d epend ent of the transit photometry - 
the iBaines et al.l (|2007t ) investigation relied upon tran- 
sit photometry for a value of the planet-star radius, R. 



There is further the possibility that conducting such an 
observation in a wavelength-dependent sense could probe 
molecular opacity effects of the the planetary atmosphere 
though sensing apparent radius dependencies. For exam- 
ple, the recent detection o f methane as an co nstituent of 
HD189733b's atmosphere (|Swain et al.ll2008f l provides a 
tantalizing goal for this technique; in principle, a suf- 
ficiently precise measurement of this nature, using nar- 
row channels inside of the iJ-band (say, comparing 0.05 
fim-wide channels centered at 1.5, 1.66, and 1.73 re- 
spectively) could confirm this detection by detecting the 
wavelength dependence of the planetary radius. How- 
ever, this appears to require levels of closure phase pre- 
cision beyond even the best cases consider ed here - an 
examination of the absorption depth data of'Swai n et al.l 
(2008) indicates such a detection to require radius mea- 
surement precision at the 0.2 — 1.1% level. 

These simulations use only the time-tagged closure 
phase data from a single transit event; supporting photo- 
metric and radial velocity signatures of the transit (and 
interferometer visibility measures), or multiple transits, 
have the potential to significantly improve the quality 
of the {r star, R,ci,b, JDq} fit parameters. Investigators 
wishing to mesh such data sets will of course have to 
pay particular attention to uniform time-tagging. It does 
seem possible, however, that data sets of such richness 
will be able to probe other system parameters: possi- 
ble moons of the transiting planet, and the presence of 
other stellar planetary companions due to variations in 
the transit timing and impact parameter. 

This approach is the only currently available technique 
that provide any value for the transit event orientation 
angle, a; it is also an independent check on parameters 
such as stellar radius or planet-star radius ratio that is 
derived from other techniques, such as spectroscopy or 
photometric timing. Use of results from such intcrfero- 
mctric observations could be highly useful for planning 
observations of TPF-I, Darwin, and other instruments 
that have a position-angle dependent response. For ex- 
ample, the extreme adaptive optics systems that are en- 
visioned carrying out planet searches an d/or character- 
izations through 'dark hole' te chniques (ISerabvn et al.l 
l2Q07t lOppenheimer et aD 120081 ) could have their search 
times reduced through a priori knowledge of the planet's 
orbital plane position angle. 
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TABLE 3 

Results from fitting transit event parameters to three synthetic transit closure phase data sets as discussed in [jg] Original values were JDo = —6000 sec, R = 0.17, 

b = 121.0/iAS, rstar = 185/iAS, and V = 0.07000/iAS/SEC. 



Input 




Time 




Ratio 


Impact 


Recovered 


Star 




Planet 


Planet 


Planet Radius 


Position Angle 




Zero 




Flanet-btar 


Parameter 


Position Angle 


Radius 




Speed 


Radius 


Fractional 


(a) [degl 


[deg] 


(Jl>0) [sec] 




Haaius [H) 


(b) [fias\ 


(a) [degj 






[v) [/ias/sec] 


\ r. 1 

(J^planet) l/^asj 


Error 


28 


0.005 


-5981 ± 20 


0. 


1709 ± 0.0061 


119.9 ± 2.6 


27.87 ± 0.26 


185.2 ± 1.7 


0, 


.07000 ± 0.00047 


31.7 ± 1.17 


27.1 


28 


0.010 


-6075 ± 38 


0. 


1701 ± 0.0074 


123.2 ± 2.2 


28.25 ± 0.42 


185.5 ± 2.7 


0, 


.06900 ± 0.00092 


31.6 ± 1.45 


21.8 


28 


0.020 


-5969 ± 53 


0. 


1687 ± 0.0115 


114.1 ± 5.9 


27.02 ± 0.82 


183.2 ± 4.3 


0, 


.06900 ± 0.00180 


30.9 ± 2.23 


13.9 


28 


0.050 


-6268 ± 97 


0. 


1717 ± 0.0126 


137.1 ± 6.3 


27.67 ± 1.40 


189.1 ± 7.3 


0, 


.07100 ± 0.00352 


32.5 ± 2.69 


12.1 


28 


0.100 


-5917 ± 183 


0. 


1748 ± 0.0352 


124.5 ± 12.2 


25.98 ± 2.98 


183.3 ± 18.2 


0, 


.07000 ± 0.00640 


32.0 ±7.19 


4.5 


115 


0.005 


-5955 ± 28 


0. 


1719 ± 0.0047 


123.1 ± 1.2 


115.52 ±0.37 


184.3 ± 1.9 


0, 


.07000 ± 0.00097 


31.7 ±0.93 


34.2 


115 


0.010 


-6020 ± 21 


0. 


1690 ± 0.0038 


120.3 ± 1.8 


114.95 ±0.42 


185.3 ± 1.5 


0, 


.07000 ± 0.00095 


31.3 ±0.75 


41.8 


115 


0.020 


-6011 ±41 


0. 


1700 ± 0.0075 


120.1 ± 3.4 


114.84 ±0.85 


184.4 ±2.4 


0, 


.06900 ± 0.00143 


31.3 ± 1.44 


21.7 


115 


0.050 


-5874 ± 77 


0. 


1708 ±0.0121 


123.9 ± 5.6 


115.37 ± 1.70 


180.3 ±5.7 


0, 


.06800 ± 0.00321 


30.8 ±2.39 


12.9 


115 


0.100 


-5957 ± 155 


0. 


1742 ± 0.0178 


136.2 ± 12.8 


115.34 ±2.64 


184.7 ±8.1 


0, 


.06900 ± 0.00507 


32.2 ±3.58 


9.0 


170 


0.005 


-6005 ± 12 


0. 


1694 ± 0.0051 


121.4 ± 1.6 


170.13 ±0.31 


185.6 ± 1.8 


0, 


.07100 ±0.00121 


31.4 ±0.99 


31.6 


170 


0.010 


-6001 ± 17 


0. 


1703 ± 0.0053 


121.8 ± 2.1 


170.05 ± 0.44 


183.8 ± 1.5 


0, 


.06900 ± 0.00065 


31.3 ± 1.01 


31.1 


170 


0.020 


-5920 ± 36 


0. 


1698 ± 0.0076 


120.2 ±4.5 


168.51 ±0.67 


181.0 ±2.9 


0, 


.06800 ± 0.00182 


30.7 ± 1.46 


21.0 


170 


0.050 


-5943 ± 65 


0. 


1715 ±0.0193 


122.0 ±6.6 


168.04 ± 1.31 


186.7 ±5.3 


0, 


.06900 ± 0.00281 


32.0 ±3.72 


8.6 


170 


0.100 


-5893 ± 236 


0. 


1696 ± 0.0275 


123.7 ±13.6 


168.12 ±4.87 


178.1 ± 13.3 


0, 


.06400 ± 0.00756 


30.2 ± 5.39 


5.6 
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